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Abstract 

The creation of particles with large transverse momenta in high energy hadronic collisions is 
a long standing problem. The transition from small- (soft) to hard- parton scattering 'high-px' 
events is rather smooth. In this paper we apply the non-extensive statistical framework to calculate 
transverse momentum distributions of long lived hadrons created at energies from low [^/s ~ 10 
GeV) to the highest energies available in collider experiments [y/s ^ 2000 GeV). Satisfactory 
agreement with the experimental data is achieved. The systematic increase of the non-extensivity 
parameter with energy found can be understood as phenomenological evidence for the increased 
role of long range correlations in the hadronization process. 

Predictions concerning the rise of average transverse momenta up to the highest cosmic ray 
energies are also given and discussed. 
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INTRODUCTION 



In what follows we write the momentum, p, in units of energy (e.g. GeV); strictly, of 
course it should be in units of energy divided by c. 

For soft (low p± events) particle production, the QCD-based string model predicts that 
the transverse momentum (transverse mass: m± = \/p±'^ + rri^) distribution of the produced 
quarks should be, in general, of the form 



d(T 



~ e-"™^'/"' , (1) 



dp_L^ 

where k is the string tension. According to Q|, when the string tension may fluctuate, the 
string hadronization becomes consistent with thermal behaviour, and 

^ ~ e-™-/^. (2) 

This means that in the soft (small p^_) region, the partonic string fragmentation picture, 
which is based on first principles, can be, to some extent, successfully replaced by the 
"phenomenological" statistical model. 

Much more interesting, and complex, is the case of high transverse momentum physics. 

The main reason for the failure of the traditional thermodynamical models of multi- 
particle production [3], was the experimentally observed significant increase of the 
production of particles with high transverse momenta discovered in the mid-seventies in 
the Intersecting Storage Ring (ISR) experiments. Its interpretation in the framework of jet 
models strongly supported the parton idea. But then new problems arose. The "hard" 
parton scattering expected from field theories implies roughly that the momenta should fall 
off as Pi.""^ 

~ F^(x„, ql) Fs{x,, ql) (3) 
dp±^ ql 

(where the F's are the respective structure functions). However the power- law fit to the 
data is closer to a fact that was realized already in, e.g., 0] and p. 

The commonly accepted explanation of high p± behaviour in inclusive pp (and pp) spectra 
is that at high enough energies the QCD hard scattering effects start to play a significant role. 
The transition "soft"— not-so-hard"— hard" is believed to combine the exponential low- 
p± domain with the asymptotic p_C'^ tail in such a way that the expected distributions are, 
around a few GeV/c and for the interaction energies available at present, just as observed. 



To follow this idea in detail further QCD calculations are needed, and to cover the whole 
p± range some interplay of perturbative (high p± - "hard" ) and non-perturbative ( "soft" ) 
physics have to be developed in a self-consistent way. The problem is not trivial and here 
we point a way forward. 

In the present paper we will show that the problem can be quite satisfactorily solved by 
way of statistical language as an effect of multi-particle, long-range correlations. We believe 
that both descriptions are equivalent. 

In the standard statistical (thermodynamical) models of hadronization, high p^'s can be 
explained (?) only by introducing additional mechanisms. Two such attempts have appeared 
recently. ^ 

One way ((9]) is to smear the pj_s with an additional, dynamic, term describing the 
collective motion of different parts of a "pre-hadronizing" state of matter. It can be achieved 
by introducing a distribution in momentum space of many fireballs created in the collision. 
However this mechanism cannot be responsible for very high p± tails. 

The second idea (Q, is to allow the Hagedorn temperature to fluctuate around 



its mean. A very good description of the data can be obtained assuming that 



with parameters 



/(;3)=p(l/r) = jgj(i)"'exp(-|) (4) 



.T/ a ' \T2/ \T/ 
It is clear that the relatively narrow F distribution in /3 = 1/T gives necessary very substan- 
tial tails in the distribution of T. The situation is presented in Fig. ^ In fact, interesting 
cases with high p_i_ originate in events with actual Hagedorn temperatures much greater than 
the "hadronic soup boiling temperature" [isj. This is rather inconsistent with the general 
Hagedorn fireball picture. 

In the present work we modify the classical Hagedorn idea by introducing the long-range 
correlations in statistical way. 

A STANDARD STATISTICAL MODEL 

One way to describe the statistical properties of the system is by introducing the concept 
of the partition function. Its classical definition for the system in the state described by the 
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P=l/T (MeV"') 



T (MeV) 



FIG. 1: Distribution of /3 = 1/T (left) and T (right) in the picture p± - broadening by temper- 
ature fluctuations with A = 10 and Tq = (l/T)"^ = a/A = 110 MeV. 



vector Qq {canonical statistical ensemble) is given by 

^(Q) = Ee-^/^^Q,Qo, 



(6) 



states 



where E is the energy of the system. For the purposes of the present paper, where the 
system contains hadrons created in high energy hadronic collisions, Q can be limited to a 
three dimensional vector Q and its components are the total charge, the baryon number 
and the strangeness of the system {Q, B, S). Of course, when needed, a generalization is 
straightforward. 

We follow here the formalism developed by Becattini (see, e.g., jl^): 

By using the integral representation of the Kronecker S factor, Eq.(jn} becomes 

2n 2n 2tt 



2{Q) = E 



states 



(2vr)^ 



(7) 





If we have a quantum gas with A*";, boson and Nf fermion species a summation over states 
can be performed, giving 



2(Q) 



1 



(2vr)^ 



^^log l-e-^-^'^^^ 

Nt k ^ 



+ EElog fl + e-^ 



Nf k 



(8) 



The sum over phase space cells, k, can be, in the continuous limit, replaced by an integration 
over momentum space: 



E 

k 



(2J. + 1) 



V 



{2ny- 



(9) 



The average multiplicity of the i^^ hadron type can be obtained from the partition function 
by introducing the fictitious fugacity factor Aj 



(n.) = ^log(Z(Qo, A.)) 



Ai=l 



thus 



{n,) = (2J, + 1) 



V 



2tt 2tt 2n 



(10) 



(11) 







where the upper sign is for fermions and the lower is for bosons. Because the e~^^'^ fac- 
tor is expected to be small for all particles except pions (T ^ 100 MeV) the following 
approximation can be made in such cases 

1 



(12) 



and then 



2n 2n 2tt 



ZiQo) i2nf 



d='</)Z(Qo) e-*'i'<^ (2J,, + 1) 



V 



ODD 



(2vr) 



d^p e 



-E/T 



^(Qo-q.) (2j^ + i)_Z_ Id'pe-^/T 



(13) 



Very good agreement with the measured particle ratios was found in, e.g., jl^. Eqs. ffTTl IT^ 
can also be used to calculate the respective transverse momentum distributions for pions 
and heavier hadrons produced in the thermodynamica]_^hadronization process. 
The high tails of these distributions are known ^] to fall like 



— — — ^ p^'/'e-^-'^ 



(14) 



what is in agreement with the low energy experimental results but there is a clear underes- 
timate of the high p_l's at collider energies. 



MODIFICATIONS OF THE STATISTICAL HADRONIZATION MODEL 

The statistical way is "phenomenological" in a sense. The conventional Boltzmann-Gibbs 
description shown above should be modified, and a new parameter describing the correlation 
"strength" , however defined, ought to be introduced. Of course, in the limit of the absence 
of correlations the new description should approach the Boltzmann form. 
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From the theoretical point of view there could be infinitely many "generalized" statistics. 
Following 'Ockham razor' we should choose that which is simple and has a clear theoretical 
background. In the present paper we test the possibility, proposed by Tsallis jl^, based on 
the modification of the classical entropy definition 

w 

Sbg = - pilnpi (15) 

i 

of the form 

w 

1 - E 

S, = (16) 

introducing the new parameter g, the non-extensivity parameter. This modification has 
been adopted in other physical applications (see, e.g., Iq). 

Maximization of the entropy requirement with the total energy constraint 

= Eo (17) 



i 

leads to the probability given by 

= ^ [1 - (1 - - E,)f^'-'^^ , (18) 

where Zq is the normalization constant related to Z{q) of Eq.® and the Boltzmann terms 
are replaced by the probabilities of the form given in Eq. ()18|) . 

It can be mentioned here, that the name of this generalization: "the non-extensive statis- 
tics" comes from the fact that the entropy Sq defined by Eq. (fTT)|) . in opposition to the Sbg, 
is non-extensive parameter (the entropy of the system consisting of two separated parts is 
not a sum of their entropies). A similar statement is valid for the total energy of the system. 
If it consists of n isolated parts (e.g., gas of "non-interacting" hadrons), each of energy Ei, 
the system total energy is not equal to Ei but given by 

E = j^^i + {q-^)/TYlE^E, + {q-lf/T'YlE^EjEk + ... . (19) 

i i,j i,j,k 

Additional terms can be interpreted as an effect of the long-range correlation which appears 
to be not only the mathematical construction adopted here, but has something to do with 
the real cause of the correlation. 

Equation (fTHj) can be rewritten introducing new symbol, e^, defined as 

(e,)^ = [1 + (1 - q)xY^^'-'^ (20) 



(for completeness, in the q=l limit we have, as we should, ef = e^). Then 

Pf. - [1 - (1 - <l)/UE, - Eo)f^'-'^ ~ ef^/^^ (21) 
and the partition function can be written in the form 

Z,(Q) = E -^Tjld^ef'^S'^^-''^^ . (22) 

states l^TTj ^ ^ ^ 

With such a modification, the transverse momentum distribution becomes 
for particles other than pions, while for pions the modified Eq.(|TT|) should be used: 

27r 27r 2tt 

I I Q Q 

(the vector in the above equations represents the i^^ particle type and q is the non- 
extensivity parameter) . 

In the present work we have evaluated Zg functions for a variety of its parameters T, V, 
and q, and for Q values which cover the production of over 100 hadrons of masses below 
2 GeV. The decays of short-lived particles were then performed (with three-body decay 
products distributed uniformly in the Dahtz plot). The effect on the p± distribution caused 
by the decays is shown in Fig. |21 

Examples of final transverse momentum distributions are shown in Fig. 01 for different 
thermodynamical parameters: T, V, and q. Closer inspection of the pj_ distributions gives 
reasons to reduce the number of (independent) parameters which should be used when 
comparing model calculations with the measured data in the high p± region. The dependence 
on the hadronization volume V for not very small p±, as seen in the left graph in Fig. El is 
not crucial. Additionally, the parameter T (hereafter it will be called 'temperature', bearing 
in mind that in the non-extensive thermodynamics, for q > 1, its meaning is not so obvious) 
and V are strongly correlated. In Ref. jl^, for the fitting procedure, the parameter VT^ has 
been chosen instead of V. The main subject of the present paper is not the multiplicity, but 
rather the shape of the transverse momentum distribution (and mainly for high transverse 
momenta), thus the normalization in the p± region of interest is quite natural to be used. The 
temperature influences mainly relatively small momenta, not much higher than its actual 



10 



04H 

Oh 



10 



10 



q= 1.05 
V=20 fm' 
T=130MeV 



10 



(GeV) 



FIG. 2: Distributions oi p± for charged primary produced hadrons (dashed hue) and for long-hved 
hadrons after all decays (solid line) (r=130 MeV, V=20 fm^ and (7=1.05). (Strictly, the abscissa, 
p±, should be in units of GeV/c but we disregard the 'c' here and elsewhere.) 

value, if, of course, the normalization is treated separately. The general normalization factor 
combines both T and V dependencies and, as detailed calculations confirm, when analyzing 
high p± data (about and above ~1 GeV/c), the parameter V can be successfully replaced by 
the overall normalization parameter with only a slight change in two remaining parameters: 
T and q. 

The non-extensivity parameter q, the crucial one for the present work, determines the 
asymptotic index of the high pj_ distribution tails. The following approximate formula found 
in Ifl illustrates this 



f{p±) ~ P± y dpi, [l + (l-g)/T 





+ P± + ^ 



{pjT f [l+p^(g-l)/T]- 



9-1 ' 2 



(25) 



We will compare it with our exact calculation results. 

Investigation of the energy dependence of q should give us the answer to the question as 
to whether the modified thermodynamical model can be reasonably applied as a phenomeno- 
logical description of complex QCD principles "at work" in the hadronization process. 
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FIG. 3: Distributions of p± for charged long-lived hadrons calculated for different V, T, and 
values around r=130 MeV, V=20 fm^ and g=1.05. 

APPLICATION TO THE DATA 



The non-extensive framework was successfully applied to transverse momentum data for 



hadrons in 



17l |. In the present work we present results concerning pp (and pp) 



reactions. The available data covered quite a wide range of particle interaction energies. We 
have started our analysis at piab = 100 GeV/c 3] where the p± distributions match quite 
nicel y th e exponential behaviour and then through [l^ and ISR energies 0, El, SPS 



22 



and [2j| up to two Tevatron energies: = 630 and 1800 GeV j2J, |25[ . 

Our model parameters (temperature T and non-extensivity q) have been adjusted (for 
p± > 0.5 GeV) to describe measured invariant p± distributions. The absence of any system- 
atic change of temperature was found. Thus we set the value of T equal to 130 MeV and 
performed the minimization procedures again. 

The first important point we have to mention here is that the reproduction of the data 
is very good. 

The next finding is a systematic increase of the non-extensivity parameter starting from 
the value of 1 (i.e., exponential, e.g., classical, p± distribution tails) at y/s ^ 10 GeV. 

The near-perfect match as obtained suggests a check to see if there exists any evidence at 
all, that the non-extensivity parameter may have a non-unique value for a fixed interaction 
energy. Some interesting statements, sugges ting such q behaviour (but for the longitudinal 



phase space) have been published recently 



26[. For example, it is possible that the non- 



extensivity varies with the actual multiplicity of created particles, if so it could be related 
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somehow to, for example, the impact parameter of the colhding hadrons |27l |. 

The spread of q was allowed to be Gaussian (with the mean value and dispersion as 
free parameters to be adjusted). Higher degree polynomials have also been tested, but no 
improvement was found. 

The fitting procedure was repeated again for all data sets listed above with the T param- 
eter released free again at the first step. No important energy dependence was noticed here 
either, so we fixed it again at 130 MeV. 

The final transverse momentum distributions are presented in Fig. |3] 

The distributions of the q parameter starts from 6{q — 1) for sub-ISR energies, and 
then, when the average q differs from 1 they become Gaussian (truncated below 1). Some 
examples are presented in Fig. El In all cases the widths of the q distributions are very small 
(their dispersions are about 0.01). This allows us to claim that the transverse momentum 
distributions for a given interaction energy can be well described by a single value for the 
non-extensivity parameter. This means that the parameter does not depend on the actual 
multiplicity (impact parameter) or any other interaction characteristic which may fluctuate 
from event to event. Its value is determined only by the available center of mass energy. It 
should be mentioned here, that this statement supports the non-extensive thermodynamical 
treatment of the hadronization process making it simple and clear. 

However, most important is the intriguing regularity as seen in the Fig. IHl where we show 
how the non-extensivity parameter q depends on the interaction energy. 

The value of q should not exceed 1.25 [l^, and it is quite reasonable to expect this as an 
asymptotic limit leading to QCD-inspired p±^ distribution for large transverse momenta at 
extremely high energies. The simple dependence is found to be of the form 

q = 1.25 - 0.33 (26) 

as shown in Fig. IHl 

RISE OF THE AVERAGE TRANSVERSE MOMENTA 

The average transverse momentum of particles created in high energy hadronic inter- 
actions in the non-extensive thermodynamical model increases when the non-extensivity 
parameter increases. In |16::] the average of the distribution given by Eq. ()25|) was calculated. 
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FIG. 4: Transverse momentum distributions measured at different energies compared witfi non- 
extensive fits. 
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FIG. 4: (cont.) Transverse momentum distributions measured at different energies compared with 
non-extensive fits. 
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FIG. 6: Energy dependence of the average value of the non-extensivity parameter q obtained from 
the p± distributions. The functional fit is shown by the line. Uncertainties at each point coming 
from the fitting procedure are much smaller than plotted points. The spread with respect to the 
line is a result of differing systematic uncertainties in the p± distributions measured by different 
experiments. 

The form of the dependence is 



This is, however, an approximate result (as is Eg. ()25|) ). The approximation is quite good 
but it can only give the average value for primary created hadrons. We have obtained the 
average p± exactly and corrected for all the effects of decays of short-lived resonances. It 
is presented in Fig. [7| where Eq. ()27|) is compared with the results obtained from the fits 
presented in Fig. HI 

Using Eq. we can present the average p± as a function of energy; the result is shown 
in Fig. IHl 

The existence of the 1.25 asymptotic limit for q and the clear and well defined rise of q 
observed in the data allows us to extrapolate the results to much higher energies with a high 
degree of confidence. The line plotted in Fig. |H1 shows a moderate {p±) increase (as s0.037 
or of about 0.1 GeV per decade of the center of mass available energy) for large interaction 
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FIG. 7: The calculated average transverse momentum of long-lived charged hadrons as a function 
of the non-extensivity parameter q. The dashed line shows the approximation given by Ea. (|27() . 
The points represent the mean p±^s calculated from the fits shown in Fig. [l] 
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FIG. 8: The average transverse momentum of long-lived charged hadrons as a function of the 
available interaction center of mass energy for our model with q{s) given by Eq. (|26|) . The points 
are as in Fig. [7| 
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energies. 

This is relevant for the creation of new Monte Carlo generators and interaction models 
used in cosmic ray physics simulations of extensive air showers (EAS) at very high energies. 
Just now is the moment when the new generation of giant EAS arrays are built and new 
data from, e.g., Auger Observatory for energies as high as 10^'^ eV are expected |2S,]. Their 
proper interpretation and the accurate primary particle energy estimation requires detailed 
knowledge of the expected lateral particle distribution which is determined in part by the 
transverse particle spread high in the atmosphere where the particles interact. The energies 
are, of course, much above those energy created at the accelerators. 

According to [29j the average p± used in all models implemented in the CORSIKA code, 
one of the most widely used EAS simulation programs, at lO^^eV is 0.55 GeV while our fit 
suggests a value which is about 20% higher. 

Cosmic ray experiments have been delivering for a very long time data related to average 
transverse momenta at energies exceeding the contemporary accelerator abilities. The most 
straightforward, at first sight, are the calorimeter experiments at mountain altitudes. The 
claim of an abrupt and very substantial p± increase at around y/s = 1000 GeV made in 2^ 
was based on observations of events with large values of the Er product in the Tien-Shan 
hadronic calorimeter. The registration of the shower of particles and determination of their 
energies E and distance to the shower axis r enables a determination of the distribution of 
product Er which is related to the transverse momentum by 

r E 

Px = ^-, (28) 

where h is the actual hadron pro duction height. However, methodological difficulties in the 
data interpretation (see, e.g, ^l|) do not permit such a strong statement. In spite of this 
fact Er is still the simplest and promising way to study the transverse spread of hadrons 
created in very high energy collisions. 



In 



32l ] the distribution of Er from the joint Chacaltaya-Pamir experiment was published. 



again with the conclusion that its shape could not be fully explained by the simulations. 
Their data are shown in Fig. 121 as solid circles and the results of our calculations with two 
interaction models are given by the respective histograms. For the shower simulation we 
used the structure of the CORSIKA program The well known FRITIOF model 
was implemented to see if the event generator widely used in accelerator physics works well 
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FIG. 9: Distribution of {Er) in the Chacaltaya-Pamir experiment [33] compared with predictions 
of FRITIOF (dashed histogram) and modified HDPM (sohd one) models 



also in the cosmic ray domain. The second model was the version of the default CORSIKA 
model called HDPM (js^ with further improvements) in which we change the generation of 
transverse momenta according to the results obtained in the present paper. 

It is seen that the data, at least at high {Er), can be explained by the simulations. The 
FRITIOF generator does not look very well here, but it do not have to be only the problem 
of transverse momenta. For the emulsion chamber data the very forward region of particle 
creation is essential and the FRITIOF was rather tuned (and the ARIADNE |3 

which is 

the part of it responsible for hard gluon bremsstrahlung) for other sorts of data. This is 
seen in FiglTUl where we plot the distribution of the energy fraction carried by the energetic 
photons (7 quanta) in the 7-hadron families observed by the Chacaltaya experiment. 



CONCLUSIONS 



We have applied modified, non-extensive statistical thermodynamics to the data on par- 
ticles with high transverse momenta created in high energy hadronic collisions. It has been 
found that the data can be successfully described with a fixed temperature parameter T and 
a non-extensivity parameter q rising slowly with interaction energy. 
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FIG. 10: Distribution of electromagnetic energy fraction seen in hadron families at mountain 







altitude cosmic ray Chacaltaya experiment |33|- The dashed histogram represents FRITIOF model 
prediction while the solid one is the modified version of HDPM algorithm. 

The non-extensive thermodynamics can be used to describe phenomenologically the 
long-range correlations (e.g., jet phenomena) without introducing particular physics. The 
achieved agreement for one q value depending only on the interaction energy confirms such 
interpretation of the non-extensivity. 

The rise of q leads to systematic rise of the average transverse momentum. The regularity 
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allows us to predict with some confidence the rise of and extrapolate it to the energies 
exceeding these available at present accelerators up to the highest cosmic ray energies. The 
importance of the exactness of such an extrapolation is clear in view of the fact that extensive 
data from the Auger Observatory are expected soon. The interpretation and understanding 
of EAS features at primary particle energies exceeding 10^°eV is a clue to the solution of 
one of the most exciting astrophysical problems. 
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